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With regard to the toroid contributions, a modified system of equations of electrodynamics mov- 
ing continuous media has been obtained. Alternative formalisms to introduce the toroid moment 
contributions in the equations of electromagnetism has been worked out. The two four-potential 
formalism has been developed for the electromagnetic continous media subjected to Lorentz trans- 
formations. 
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^ PACS 03.50.De Maxwell theory: general mathematical aspects 
^ . PACS ll.lO.Ef Lagrangian and Hamiltonian approach 



1. INTRODUCTION 



The history of electromagnetism is the history of struggle of different rival concepts from the very early days of 
', its existence. Though, after the historical observation by Hertz, all main investigations in electromagnetism were 
' based on Maxwell equations, nevertheless this theory still suffers from some shortcomings inherent to its predecessors. 
Several attempts were made to remove the internal inconsistencies of the theory. To be short wc refer to very few of 
them. One of the attempts to modify the theory of electromanetism was connected with the introduction of magnetic 
charge in Maxwell equation by Dirac while keeping the usual definition of E and B in terms of the gauge 

potentials. A very interesting work in this direction was done by Miller [Q. In this paper he showed the mutual 
substitution of the sources as follows: 



, ^ p^w^-f ^ ccurlM ^ -divM = ^ p"V = j™ = ccurlP -> -divP = p" ^ 

^ Recently D. Singleton |^,^ gave an alternative formulation of classical electromagnetism with magnetic and electric 
O charges by introducing two four- vector potentials A^^ — ((/)e. A) and — (0m, C) and defining E and B fields as 



E = -V0e - ^ -curie (1.1a) 

dt ^ ' 

dC 

B = - ^ + curl A (1.1b) 

Inserting these newly defined vector potentials into the generalized Maxwell equations |^ 

9E 

curlB = — +Je (1.2a) 

dt ^ ^ 

divE = pe (l-2b) 

-curlE = — +J„ (1.2c) 
ot 

divB = p,„ (1.2d) 
where pm and are the magnetic charge and current, respectively, and imposing the Lorentz gauge condition 

^-|-divA = 0. ^+divC = (1.3) 
ot ot 

he arrived at the wave form of Maxwell's equations with electric and magnetic charges. Note that, a similar theory 
(two potential formalism) was developed by us few years ago (we will come back to it in Sec. 3). Though this 
treatment avoids the use of singular, nonlocal variables in electrodynamics with magnetic charge and makes the 
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Maxwell system more symmetric, since both the charges in this approach are gauge charges, the main defect of this 
theory in our view is that the existence of magnetic charge still lack of experimental support, hence can be considered 
as a mathematically convenient one only. 

Recently Chubykalo a.o. made an effort to modify the electromagnetic theory by invoking both the transverse and 
longitudinal (explicitly time independent) fields simultaneously, thus giving an equal footing to both the Maxwell- 
Hertz and Maxwell-Lorentz equations To remove all ambiguities related to the applications of Maxwell's 
displacement current they substituted all partial derivatives in Maxwell-Lorenz equations by the total ones, such that 

l^l-(V.grad) (1.4) 

where V = dYq/dt and Yq{t) are the velocity and the coordinate of the charge, respectively, at the instant t. Further 
they separated all field quantities into two independent classes with explicit {}* and implicit {}o time dependence, 
respectively. Thus, the component Eq of the total electric field E in every point is understood to depend only on the 
position of source at a given instant. In othe r wo rds, Eq is rigidly linked with the location of the charge. From this 



point of view, the partial time derivative in (1.4) must be related only with the explicit time-dependent component 



E* whereas the convection derivative only with Eq 

rfE dE 



dt dt 



(V.grad)Eo, E = E*(r, i) Eo(R(i)) (1.5) 



where r is a fixed distance from the origin of the reference system at rest to the point of observation and R(t) = r—rq{t). 

Another attempt to modify the equations of eletromagnetism is connected with the existence of the third family of 
multipole moments, namely the toroid one. This theory was developed by us during the recent years. Recently we 



introduced toroid moments in Maxwell equations exploiting Lagrangian formalism |10|. In the Sec. 2 of this paper 
we give a brief description of this formalism. Moreover, here we develop an alternative method to introduce toroid 
moments in the equation of electromagnetism. In Sec. 3 we develop two potential formalism suggested by us earlier. 



2. INTRODUCTION OF TOROID MOMENTS IN THE EQUATIONS OF ELECTROMAGNETISM 



In early fifties, while solving the problem of the multipole radiation of a spatially bounded source, Franz and 
Wallace 0,11^ found a contribution to the electric part of radiation at the expense of magnetization. Further 
Ya. Zel'dovich |13| pointed out the non-correspondence between the existence of two known multipole sets. Coulomb 
and magnetic, and the number of form-factors for a spin - ^ charged particles. Following the parity non-conservation 
law in weak interactions Zel'dovich suggested a third form-factor in the parametrization of the Dirac spinor particle 
current. As a classical counterpart of this form- factor he introduced anapole in connection with the global electromag- 
netic properties of a toroid coil that are impossible to describe within the charge or magnetic dipole moments in spite 
of explicit axial symmetry of the toroid coil. In 1965 Shirokov and Cheshkov |14| constructed the parametrization for 
relativistic matrix elements of currents of charged and spinning particle, which contain the third set of form-factors. 
Finally, in 1974 Dubovik and Cheskov [|l5| determined the toroid moment in the framework of classical electrodynam- 
ics. Note that anapole and toroid dipole are not the different names of one and the same thing. They are indeed quite 
differnt in nature. For example, the anapole cannot radiate at all while the toroid coil and its point-like model, toroid 
dipole, can. The matter is that the anapole is some composition of electric dipole and actual toroid dipole giving 
destructive interference of their radiation. Recently a principally new type of magnetism known as aromagnetism 
was observed in a class of organic substances, suspended either in water or in other liquids p^ . Later, it was shown 
that this phenomena of aromagnetism cannot be explained in a standard way, e.g., by ferromegnetism, since the 
organic molecules do not possess magnetic moments of either orbital or spin origin. It was also shown that the origin 
of aromagnetism is the interaction of vortex electric field induced by alternative magnetic one with the axial toroid 
moments in aromatic substances ]l7t . 

In a recent work Dubovik and Kuznetsov ||l8|| calculated the toroid moment of Majorana neutrino. It was also 
pointed out that the magnitude of the toroid dipole moment of a Dirac neutrino (j/^) is just the half of that of a 
majorana one (fm) and both of them posesses non-trivial torid moments even if = [ p^ . 

The latest theoretical and experimental development force the introduction of toroid moments in the framework 
of conventional classical electrodynamics that in its part inevitably leads to the modification of the equations of 
electromagnetism and the equations of motion of particles in external electromagnetic field. In the two following 
subsection we give two alternative schemes of introduction of toroid polarizations in the electromagnetic equations. 
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To begin with we write the Maxwell equations for electromagnetic fields in vacuum, in the presence of extraneous 
electric charge p and electric current, that is, charge - in - motion, of density j. 

laE 47r. ,^ ^ 

curlB ^ = — j 2.1a 

c at c 

divE = 47rp (2.1b) 
1 

curlE+-— =0 (2.1c) 
c at 

divB = (2. Id) 

where E and B are the flux densities of electric field and magnetic induction, respectively. Note that, the electric 
charges and electric currents, being destributed in vacuum, construct the electromagnetic structure of matter [ pO[ and 
are relateted to the elementary charge e in the following way 

p(r)=^e„(5(r-q„) (2.2a) 

n 

j(r)=^e„q„'5(r-qj (2.2b) 



So, to describe the system as a whole the equations (2A) and (2^) should be supUemented by the equation of motion 
of micro-particle, i.e., 

g 

mnq„ = e„E + — q„ x B (2.3) 

c 

In this process there occurs a vast number of problems connected with the different areas of this vast field: 

Write the classical and quantum analogies of the equations of motion of a point-like particle possessing toroid dipole 

(with usual properties); 

Solve the boundary- value problems for the model with all the of vector polarizations; 

Formulate the electrodynamics of continuous infinite media for the latter case. 

Sin ce th e electromagnetic field in media generates bound charges and bound currents, the source parts in ( p. la ) 



and (2.1b) should be supplemented as follows pi[ 



Ptotai = P + Pbound = P - divP (2.4a) 
Jtotai = J = j + ^ + ccurlM (2.4b) 

Then we may write the Maxwell equations for electromagnetic fields in media as 

,„ laD 47r. , , 

curlH--— = — j 2.5a 

c at c 

divD = 4ttp (2.5b) 
1 (9B 

curlE+-— = (2.5c) 
c ot 

divB = (2.5d) 

where we denote D = E + 4ttP and H = B — Att'M. arc the flux densities of electric displacement vector and magnetic 
field, respectively. 

Note that dP/dt, known as bound-charge current density is independent of the details of the model and is a 
conduction current. The difference between an "ordinary" conduction current density and the current density dP/dt 
is that the first involves free charge p (foreign charge over which wc have some control, i.e., charge that can be added 
to or remove from an object) in motion, when the second hound charge (the integral parts of atoms or molecules of the 
dielectric) in motion. Another obvious practial distinction is that it is impossible to get a steady bound charge current 
that goes forever unchanged. The hound currents ccurlM are associated with molecular or atomic magnetic moments, 
including the intrinsic magnetic moment of particles with spin, whereas free currents j are "ordinary" conduction 
currents flowing on macroscopic paths and can be started or stopped with a switch and measured with an ammeter. 

The foregoing equations can be obtained from a Lagrangian describing the interacting system of electromagnetic 
field and non-relativistic test charged particle [E2[ 
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L = imq^ + -L (curlA)21dr + - ( J(r) • A(r)dr (2.6) 

2 Stt J c'' c J 

Here the total current of the system J coincides with the current j = eg, generated by the free charges (charged 
particles) in motion since there is no polarization current. 

It can be easily verified that variation of the Lagrangian ( |2.6|) wi th respect to the particle coordinates gives the 
second law of Newton with the Lorentz force, i.e., the equation (|2. 



mq = eE(q, t) H q x B(q, t) 

c 



whreas the variation with respect to field variables gives (2.1a), i.e., 



l^E 47r. 

curlB--— = j 2.7 

c ot c 



whre we denote B = curlA and E = —c^^{dA/dt). It should be emphasized that E in (2^) and (2/7) is the 
transverse part of the total electric field. The lon gitudinal electric field in question is entirely electrostatic. The 
Hamiltonian, corresponding to the Lagrangian (|2.6|) reads 



H[U,A;p,q]=pq + jUAdr-L (2.8) 

^ i [p ~ ' + ^ / + 

where the corresponding conjugate momenta are p — mq + (e/c)A(q, t), n(r) — (47rc^)^^A. 

It is well known that in classical dynamics the addition of a total time derivative to a Lagrangian leads to a new 
Lagrangian with the equations of motion unaltered. Lagrangians obtained in this manner are said to be equivalent. 
In general, the Hamiltonians following from the equivalent Lagrangians are different. Even the relationship between 
the conjugate and the kinetic momenta may be changed |^^. Let us now construct an equivalent Lagrangian to that 
of ( |2.6D in the following way 

Li=L-~ J P(r) • A(r)dr - J div(M x A)dr (2.9) 

The new Lagrangian is a function of the variables q, q and a functional of the field variables A, A, and the equations 
of motion follow from the variational principle. We have 

dL 1 

^^-(J-P-ccurlM)^=j. 

^l-W^) ^ iMB-4.curlM]^ - ^[curlH], 

d dLi 1 dT>i 

= — E + 47rP. = -- 



dtd{dA^/dt) 47rc5r 47rc dt 

Applying the Euler-Lagrange equations of motion 

d dLi d dLi dLi 

dt dA,{v) ^ ^d^jd{dA,/dxj) "dA," 



we obtain (|2.5a| ), i.e.. 



ldT> 47r. 

curlH = — — H J 

c ot c 



The corresponding Hamiltonian now can be written as 



i7i[n,A;p,g] = -^[p--A(q,t)]' + ^ /[D^+B^jdr (2.10) 
2m c sn J 

D • Pdr + 27r / P^dr - / M • Bdr 
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In the same way we can introduce the toroid polarizations in the system. If we take into account that the total current 
now has the form IM 



Jtotal 



dP ^ dcuAT" 
'dt ^ 



dt 



c curlM + c curl curlT" 



(2.11) 



we again obtain ( |2.5a| ) with D = E + 47r(P + cur IT") and H = B - 47r(M + curlT™). Here T"" and T" are the toroid 
polarizations of electric and magnetic type respectively. In terms of D and B the system (2.5) can be written as 
follows 



1 dV A-K 

curlB = + 47r{curlM + curl curl T™} H j 

c dt c 

divX> = 47r/9 

curlX> = --^ + 47r{curlP + curl curl T""} 
c dt 

divB = 



(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 



where we forced to redefine B and D as B and T> respectively. Indeed, due to the introduction of toroid polarizations, 
having independent origin in terms of atomic and molecular current and charge destributions, the quantities B and 
D as well as E and H lost their initial meaning. The existence of the vorticities T*^ and T™, generally speaking, can 
be imputed to the one and the same physical volume. So in what follows, we substitute E, B, D, H by £, S, X>, 7i. 
The equation of motion ( 2.3) should also be rewritten as follows 



77iq = e£ + -qx B 



The equation (2.12a) can be derived as earlier by constracting a Lagrangian equivalent to Li such that 

L2 = Li--— f curlT'=(r) • A(r)dr - / div(curlT'" x A)dr 
cat J J 

The corresponding Hamiltonian reads 

H2 [n, A; p, g] ^ [p - £ A(q, t)]' + ^ [ [D^ + B^] dr 

X> ■ [P + curlT'^Jdr + 27r / [P + cuAT^j^dr - /[M + curlT"] ■ Bdr 



(2.13) 



(2.14) 



(2.15) 



Note that in our previous work we impose the fohowing additional condition 



™ „ 1 • e,m 

rlT™-'' = ±-T 



(2.16) 



The relation (2.16) demands some comments. Both T"^ and T™ represent the closed isolated lines of electric and 
magnetic fields. So they have to obey the usual differential relations similar to the free Maxwell equations ^,^). 
However, remark that signs here are opposite to the corresponding one in Maxwell equations because the direction 
of electri c dip ole is accepted to be chosen opposite to its inner electric field ||27 ]. It should be emphasized that the 
relation (2.16) is a local one and it is not necessary to demand the condition (2.16|) to be held to introduce toroid 



polarizations to the electro magne tic equations. We would also like to remark that the equations ( 2.12a ) and ( 2.12c ) 
can be derived strait from (2.5a) and ( ^.5c ) making the following substitutions in them: 



P + curlT'^ 



M 



M + curlT" 



and redefining the vectors B and D as earlier. 



3. TWO POTENTIAL FORMALISM 



It is generally inferred that the divergence equations of the Maxwell system are "redundant" since they are the 
consequences of curl equations under the condition of coninuity [Q. Recently Krivsky a.o. claimed that to 

describe the free electromagnetic field it is sufficient to consider the curl-subsystem of Maxwell equations since the 
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equalities divE = and divB = are fulfilled identically. Contrary to this statement, it has been proved that the 
divergence equations are not redundant and that neglecting these equations is at the origin of spurious solutions in 
computational electromagnetics [^,^. Here we construct generalized formulation of Maxwell equations including 
both curl and divergence subsystems. In this section we develop two potential formalism (a similar formalism was 
developed by us earlier with the curl-subsystem taken into account only). Note that in the ordinary one potential 
formalism (A, (^) the second set of Maxwell equations are fulfilled identically. So that all the four Maxwell equations 
bring their contribution individually, in our view, one has to rewrite the Maxwell equation in terms of two vector and 
two scalar potentials. To this end we introduce so-called double potential |36 3^,|l^ to the system ( 2.12 ), i.e.. 



curlB = — — 4- 

c ot 

divX> — Airp 

curlX> = — 

c dt 

divB = 



47r . 

— Jfrc 



47r{curlM + curl curl T'"} 



47r{curlP + curl curl T'^} 



(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 



Before developing the two potential formalism we first rewrite system ( |2.1| ) in terms of vector and sc alar potentials 
A, (f) such that B = curlA, E — — V0 — (l/c)(9A/9i). Following any text book we can write system (2.1) as 



47r. 47rr. 5P 

° A = Jtot = jf,^^ — 

c c at 

□ (/) = -Ati[p - divP] 
under Lorentz gauge, i.e., div A + {l/c){d(j)/dt) = and 



;curlM 



□A 



47r r. 

[Jtot 



47r dt^ 



-47r[/9 - divP] 



(3.2a) 
(3.2b) 

(3.3a) 
(3.3b) 



under Coulomb gauge, i.e., divA = 0. Here □ = — {l/c^){d'^ /dt^). Note that to obtain (3.2) or (3.3) it is sufficient 
to consider (2.1a) and ( p.lb| ) only since the two others are fulfilled identically. 

Let us now develop two potential formalism. Two potential formalism was first introduced in p6[ | and further 
developed in [^,0. In both papers we introduce only two vector potentials a™, af^ and use only the curl-subsystem 
of the Maxwell equations with the additional condition diva™'^ = 0. Thus, in our view our previous version of 
two potential formalism lack of completeness. In the present paper together with the vector potentials q;'",q:'^ we 
introduce two scalar potentials (p™ and (f'^ such that 



B = curia" 



1 da" 



c dt 

X> = curia'' — - 

c dt 



(3.4a) 
(3.4b) 



It can be easily verified that system of equations (3.1) are invariant under this transformation and take the form 



An 

□ a™ = [j -f c curl M + c curl curl T™] , 

c 

□ yj™ 

An r 

Dcx" = curlP 

c 



curl curl T*"], 



under div ct' 



(l/c)(c)(^^''"/9t) = and 

ATT r 

□ a™ = j-f ccurlM- 

c 



: curl curl T" 



Att dt 



V2 







Att 



[curl P -t- curl curl T'^ 



1 av?" 

'A^~dt 



'At: p 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 

(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 
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under diva'"''^ = 0. The solutions to the systems (3.5) and (3.6) can be written as follows (see for example [ |lO| , |3q] ): 
The solutions to the d'Alembert equation 



look 



F{r,t) = - 



DF{r,t) ^ f{r,t) 
f{r',t')dr' 



in 



all space 



t'=t-|r-r'l/c 



whereas the solutions to the Poisson equation 



read 



F(r) 



1 

An 



fir) 
/(r')dr' 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



Let us rewrite the equations (|3.5[) and (3.6) for the fields subject to Lorentz transformation. If the fields in 



stationary frame (unprimed) are connected with those in moving one (primed) in the following way 



B 


= 7(B' 4 


-f3 X E') 


E 


= 7(E'- 


c/3 X B') 


P 


= 7(P' + 


-13 X M') 

c ' 


M 


= 7(M' - 


-cfBx P') 




= 7(T'='- 


f i/3 X T"') 

c 


rprn 


= 7(T"' 


- c/3 X T-^') 


a" 


- 7(«™' 






= 7(^'"' 


- /3a^') 


a" 




- /3^™') 






h /3a™') 


P 


= 7(p' + 


-/3J') 

c 


J 


-7(J' + 


c/3p') 



The equations (3.5) now can be written as 
47r72 



□ 93™' = /3na'= 



[7-ij' + ccurl (M' + curlT™') - curl (/3 x P' + curl (/3 x T^'))] , 



4-717^ 



nip"" = -/3na'^ 



[curi (P' + curlT'^') + -curi (/3 X M' + curl (/3 x T™'))] , 
An{p' + (33' /c) 



(3.11a) 
(3.11b) 
(3.11c) 
(3. lid) 
(3.11e) 
(3.11f) 
(3.11g) 
(3.11h) 
(3.11i) 
(3.11j) 

(3.11k) 
(3.111) 

(3.12a) 
(3.12b) 

(3.12c) 
(3.12d) 



It is necessary to emphasize that the potential descriptions electrotoroidic and magnetotoroidic media are completely 
separated. The properties of the magnetic and electric potentials a™ and a*^ under the temporal and spatial inversions 
are opposite ||2l. The potential a" (a™) is related to the toroidness of the medium T"" (T") as B (D) to M (P). 
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4. CONCLUSION 



The modified equations of electrodynamics has been obtained in account of toroid moment contributions. The two- 
potential formalism has been further developed for the equations obtained. Note that introduction of free magnetic 
current and magnetic charge in the equations (3.1c) and ( p. Id ) respectively leads to the equations obtained 
by Singleton 
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